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The quantum ratchet current is studied in the parameter space of the dissipative kicked rotor 
model coupled to a zero temperature quantum environment. We show that vacuum fluctuations 
blur the generic isoperiodic stable structures found in the classical case. Such structures tend to 
survive when a measure of statistical dependence between the quantum and classical currents are 
displayed in the parameter space. In addition, we show that quantum fluctuations can be used to 
overcome transport barriers in the phase space. Related quantum ratchet current activation regions 
are spotted in the parameter space. Results are discussed based on quantum, semiclassical and 
classical calculations. While the semiclassical dynamics involves vacuum fluctuations, the classical 
map is driven by thermal noise. 
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I. INTRODUCTION 

A general description of unbiased transport of particles 
in nature, which is named ratchet effect, is a challenging 
issue with implications in distinct areas, such as in molec¬ 
ular motors in biology [T] , nanosystems like graphene P] , 
control of cancer metastasis [2| , micro and nanofluids [3] , 
particles in silicon membrane pores [5], cold atoms [6j, 
solids and drops transport using the Leidenfrost effect 
00, quantum systems EHH], among many others. Cer¬ 
tainly the most relevant and common goal in describing 
these ratchet systems is to unveil how to control and at¬ 
tain an efficient transport by choosing the appropriate 
physical parameter combination, like temperature, dis¬ 
sipation, external forces etc. It is known m that in 
classical conservative systems the parameters must be 
chosen so that the underlying dynamics presents a mix¬ 
ture of regular and chaotic motion. In dissipative inertial 
systems the Classical Ratchet Current (CRC) is more ef¬ 
ficient when parameters are chosen inside the Isoperiodic 
Stable Structures (ISSs), which appear in the parameter 
space of ratchet models mm- Such ISSs are generic 
Lyapunov stable islands with dynamics globally struc¬ 
turally stable and come along in many dynamical systems 
H5HT2I- When stochastic effects are included, e. g. ther¬ 
mal fluctuations, these ISSs start to be destroyed and 
become blurred, even though they remain resistant to 
reasonable noise intensities (50]. This means that in real¬ 
istic systems the CRC is more efficient when parameters 
are chosen inside the ISSs. In m it was also shown 
that in certain multistability scenarios the current may 
actually be thermally activated. 

The natural question now is whether the same general 
statements can be made regarding the relation between 
the system’s parameters (or the ISSs) and the Quantum 
Ratchet Current (QRC). With the ongoing technologi¬ 


cal developments engineering smaller and smaller devices 
this question becomes crucial for the observation of di¬ 
rected transport in quantum systems. It has been shown 
[2l] that for specific points in the parameter space the 
QRC apparently has the same main properties of the 
CRC inside the ISSs when considered with a certain finite 
temperature (see also jB]). In this context one expects 
the quantum version of the ISSs to become blurred and 
gradually disappear as the full quantum limit is reached. 
In principle though, the simple addition of thermal ef¬ 
fects in the classical dynamics is certainly not enough to 
reproduce the quantum dynamics. 

In the present work we calculate the QRC using a 
dissipative zero temperature master equation. There¬ 
fore, fluctuations are of quantum origin while no ther¬ 
mal noise is considered in the description. In the semi¬ 
classical limit, the master equation can be recast as a 
semiclassical map, allowing us to study in details the 
quantum-classical transition. The semiclassical results 
are compared to the classical map for the ratchet system, 
derived through direct integration of the Langevin equa¬ 
tion over a kicking period. The classical map involves 
correlated thermal fluctuations. We show how and to 
which extent the quantum-classical transition affects the 
optimal currents inside the ISSs. Even though vacuum 
fluctuations destroy the well defined shape of the ISSs 
when the current is plotted, the borders of the structures 
tend to survive when the quantum and classical currents 
are statistically compared (22] . In addition, the quantum 
fluctuations can assist current activation and we spot the 
activated regions in the parameter space. 

The paper is structured in the following way. In Sec.|TT| 
we present the quantum, semiclassical and classical mod¬ 
els used to investigate the ratchet current. In Sec. |III| the 
results for the quantum ratchet current, its activation, 
and the distance correlations are shown in the parameter 
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space and discussed. Section |IV| presents the summary 
of our findings and motivates related experiments. 


II. MODELS 

A. The Quantum Problem 

The quantum problem can be described using the 
kicked Hamiltonian in dimensionless units 
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where (x,p) are position and momentum operators, m = 
1,2 ,N represents the discrete times when the “kicks” 
occur, r is the kicking period, and K is the parameter 
which controls the intensity of the kick. Dissipation is 
introduced between the kicks by coupling the ratchet sys¬ 
tem to a zero temperature environment. In usual Born- 
Markov approximation a master equation for the reduced 
dynamics is obtained [23j. Here we determine the corre¬ 
sponding density operator as an ensemble mean over pure 
states obtained from the corresponding quantum state 
diffusion [53j Ito-stochastic Schrodinger equation 

I dip) = -iH\ip)dt + ^ (L^ — (L/j,)) \i/j)d^ 

( 2 ) 

-\ e ( 4 ^ - +k^>i 2 ) 

a 


demanding since more and more states and a huge num¬ 
ber of time steps must be taken into account in this limit. 
To overcome this difficulty we use the semiclassical map 
proposed years ago |53] for the standard map. In our 
ratchet case the semiclassical map at zero temperature 
takes the form 


Pm+ 1 = IPm + 7 K { sin (dm + 4+1) 


+ - sin 


2(<?m T 'fm- f-l) T 2 f f 1) (3) 


(1 — 7 )t _ - 

Qm+1 — Qm T q i Pm +1 T tym+li 

7l»nl 


where (x m ,p m ) are scalars being the position and mo¬ 
mentum at discrete times m = 1,2,..., TV, and K is the 
kicking parameter. The map p| and the noise sources 
(■i/’m+i) Vm+i) arise from conveniently writing the propa¬ 
gator of the Wigner function as a quasi-stochastic map 
(for more details see (23j). The noise terms are uncorre¬ 
lated in time and satisfy (rj m+ i) = (ip m+ i) = 0, and 
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where (.) stands for the expectation value. The Lindblad 
operators are given by L\ = g ]>f n \Jn + 1 |u) (n + 1| , and 
L 2 = g n \/n + 11 —n) (—n — 1 |, with g being the cou¬ 
pling constant and n = 0,1,.... The | n) states are the 
eigenstates of the momentum operator p\n) = 27rn|n). 
The Lindblad operators induce a damping -—A(p), with 
rate A = g 2 . The semiclassical limit of Q and |2| is 
obtained by taking r —> 0 and A -A 00 , such that the dis¬ 
sipation parameter 7 = e~ Ar remains constant. It turns 
out that h e ff = t plays the role of the effective Planck 
constant (6 . The QRC is obtained from Cq = ||(p)||, 
where ||.|| is a double average, over quantum realizations 
and time. As initial condition (IC) we use a coherent 
wave packet localized inside a minimum of the ratchet 
potential, i. e. (£ 0 ) = 7t/ 2 with (ho) = 0.0. 

B. The Semiclassical Problem 

We can better understand the quantum-classical re¬ 
lation regarding the directed transport by studying the 
QRC as we approach the classical limit h e // -A 0. It 
is not reasonable, however, to perform a full quantum 
calculation for this purpose because it is numerically too 


Higher order cumulants were neglected (higher orders in 
h e //) so that the map ([ 3 ]) can be interpreted as a classical 
map with Gaussian noise terms of quantum mechanical 
origin. This map is valid as long h e // <C 1. The current 
obtained from pi is referred to as SRC (Semiclassical 
Ratchet Current). ICs with ((xo)) = {{po)) = 0 were 
taken inside the unit cell (—27 t, 27t). 


C. The Classical Problem 

The connection between the dynamics of a quantum 
system at zero temperature and its finite temperature 
classical limit in the parameter space is not obvious at 
first sight. Although one could argue that the first-order 
(in h) quantum correction to a classical system is equiva¬ 
lent to a white noise term |25| , the fluctuation-dissipation 
relation can differ, and it is also not clear whether a cer¬ 
tain region of the parameter space would require higher 
order corrections for a reliable description of the system’s 
dynamics. We access this connection by means of a direct 
comparison between the calculated current in parameter 
space for the quantum system and its classical finite tem¬ 
perature limit. This limit is obtained via integration of 
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the kicked Langevin equation 


p = —Xp + K 
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where r c is the classical kicking time, a = \Jl\kBT , 
and £(t) is the white noise satisfying (£(f)) = 0 and 
(£W£( S )) = — s). The integration is performed be¬ 

tween two kicks of the potential including the kick in 
the beginning, but not the one at the end. More pre¬ 
cisely, the integration is performed between to = m, — e 
and t = m + 1 — e, with positive e —> 0. Defining 
7 = exp(—A t c ), the integration produces the following 
map 
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properties (a m ) = (fi m ) = 0 , and 
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with Te// = ksTTp. In the map, the variables x m and 
p m represent respectively the position and the kinetic 
momentum of the particle when it is kicked by the m-th 
time (m = 1, 2,..., N). Notice that the correlated thermal 
noise in the map derived from the kicked Langevin 
equation with white noise, differs from the uncorrelated 
noise used in previous works 1201121 ]. As in the semiclas- 
sical case, ICs with ((x 0 )) = {(po)} = 0 were taken inside 
the unit cell (—27r, 27 t). 


III. RESULTS 


A. Ratchet current 

Now we can discuss results by comparing the cur¬ 
rents of the three problems above. Figure [T] shows the 
ratchet currents (see colors) (a) Cq (quantum), (b) Cs 
(semiclassical) and (c)-(d) Cc (classical) in the parame¬ 
ter space (A'', 7 ) with K' = (1 — 7 )tAT/| I 117 I being the 
rescaled kick parameter. In order to compare Cs and Cc 
with the Cq, we used Cs = (1 — 7)' 7 "/(7|l rL 7l)l|p|I anc ^ 
C c = (1 — 7 )r/( 7 | ln 7 |)||p||, where the double average is 
over ICs and time. We checked that the Cq converges, 


for the whole considered parameter space, after 10 3 kicks 
and some tens of stochastic realizations. For the semi¬ 
classical and classical simulations we used 10 4 iterations 
and 10 3 ICs with zero average in position and momen¬ 
tum. It is worth to mention that for the parameter 7 , 
shown in Fig.JTJa), the damping A is approximately inside 
the interval [2.0, 9.7], which represents strong dissipation. 
Thus the quantum dynamics is incoherent after the very 
short classical times m r ~ 1/(1 — 7 ) [23] and coherence 
effects are not expected to be relevant for the results 
shown in the parameter space from Fig. [lja). Agreement 
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Figure 1. (Color online) The parameter spaces showing 
the (a) Cq, (b) Cs, for h e ff = 0.082, and (c) Cc for 
T eff = 0 .082. Red (green) to yellow (white) colors are re¬ 
lated to increasing negative (positive) currents [marked with 
(—) and (+) for the printed grayscale]. Note the black bound¬ 
aries delimiting positive and negative currents. I 11 (d) we have 
the Cc for T e // = 0 with black colors related to close to 
zero currents; green to blue colors (light to dark gray) are 
related to increasing positive currents while red, yellow to 
purple (white to light gray) colors related to increasing neg¬ 
ative currents. Letters denote the chaotic background A and 
the ISSs Bi,B 2 ,C_i and F. 

between Cq [Fig. HJa)] and Cs [Fig. [TJb)] is astonish¬ 
ing, meaning that the semiclassical map ([3| nicely repro¬ 
duces the quantum results in the quantum-classical tran¬ 
sition for this range of parameters. Note that the Cc 
in Fig. |Tj(c), obtained using one fixed T e ff in the classi¬ 
cal map Q, also shows a qualitative agreement with the 
Cq. 

To understand the physical origin of the distinct cur¬ 
rents in Figs. |T](a)-(c) we analyze the Cc at zero tem¬ 
perature T e ff =0. This is shown in Fig. EL) . Two main 
regions with distinct behaviors can be observed: first, the 
“cloudy” background, identified as A in Fig. B d ) show¬ 
ing a mixture of zero, small negative and positive cur¬ 
rents, where the dynamics is chaotic M- Second, the 
four ISSs, Ai,I? 2 , C-i and F. Besides F, which has zero 
current, the ISSs are responsible for the optimal Cc and 
can be recognized by their sharp borders (for a detailed 
explanation see mm )• We note that at the crossing 
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regions of B\ and F, positive currents are observed as 
a consequence of multistable attractors. The portion of 
the parameter space shown in Fig. [T]can be regarded as 
representative of the highly dissipative ratchet regime of 
this system, presenting the main properties found in the 
larger parameter space, and is thus suitable for the pur¬ 
pose of the present work. Comparing Figs. Ba)-(c) to 
Fig- [3 d ), we can identify that the positive and negative 
C Q , C s and Cc {T e ff > 0) roughly follow the overall 
chaotic currents from the region A, but are enhanced in¬ 
side the ISS B 1 from the classical case with T e ff = 0. 
This shows that vacuum fluctuations tend to blur the 
classical ISSs, as suggested recently m- One can re- 



Figure 2. (Color online) Plotted is (a) |Cs| — |Cq| and (b) 
|Cc| — |C Q | from the data of Fig. [I] [Symbols (—)/(+) to 
identify the negative/positve currents in printed grayscale], 

cover many of the Cq’s features in parameter space using 
classical calculations. This was also observed in [23], for 
specific points in the parameter space, in a system in con¬ 
tact with a distinct thermal bath with uncorrelated noise. 
However, while a purely classical map with thermal fluc¬ 
tuations (and single T e ff value) leads to a Cc that can 
qualitatively mimic Cq in parameter space, it is usually 
quantitatively very far off. Much better agreement, both 
qualitatively and quantitatively, is reached using a semi- 
classical map. This is shown in Fig. [3 We stress that 
quantum results are very well represented by the semi- 
classical map without any free parameter. By contrast, 
the simulation in terms of the classical approach @ , in¬ 
volves T e ff as an additional free parameter that needs 
to be adjusted for every point in the parameter space to 
achieve an overall agreement with the Cq. Thus, it be¬ 
comes clear that a single T e ff is not able to reproduce 
the Cq for the whole parameter space. 


B. Correlation between SRC and CRC 

A result for the quantum-classical transition regarding 
ratchet currents is found by analyzing the correlation be¬ 
tween the Cs (at h e ff ^ 0) and the Cc (at T e ff = 0). 
It is a measure of statistical dependence, called distance 
correlation |22| . between both currents, which is zero if 
and only if the currents are statistically independent. 
Thus, the purpose of the present analysis is to search 
to which extent the quantum current is statistically sim¬ 
ilar to the classical current. For this we use the distance 


correlation, which is obtained from the expression f22] 
V(p,p) = d(p,p)/^/ A(p)A(p), where 


E SiA* (5) 
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are the distance variances. Here p refers to the time 
sequence Pi,p 2 ,p 3 , ■ ■ ■ ,Pm obtained from the map (§, 
leading to the Cs- p refers to the time sequence 
Pi,P 2 ,P 3 , ■ ■ ■ ,Pm obtained from the map Q, leading to 
the Cc- S and C are matrices defined through Sy*, = 
Sj,k — Sy — T and Cj,k = Cyfc Cy c_£, -f- c.., re¬ 
spectively. Here .syy = [p • — p k | is the Euclidean norm of 
the distance between the momenta (averaged over ICs) 
at times j and k. Thus, p ■ is an average over ICs at times 
j. The same is valid for Cy*, = |p- — p k \. In addition, Sy 
is the j- th row mean, s,k the fc-th column mean, and s,, 
is the mean value between sy and s,k (the same notation 
applies to c). Thus, each element of the matrices S and 
C contains information about: distances in time of the 
currents (already averaged over ICs) and also time aver¬ 
ages of these distances over the time series p and p. The 
distance correlation V(j),p), on the other hand, analyses 
the statistical independence of both matrices. It is not 
a direct quantity to be measured in one simulation (or 
experiment), but compares the currents obtained by two 
simulations performed in distinct semiclassical regimes. 

In Figs. [3 (a) and (b) the quantity T>(p,p) is displayed 
(colors) in the same parameter space from Fig. [3 White, 
yellow, red to blue (white, light to dark-gray) colors are 
related to increasing values of the distance correlation. 
For each point in the parameter space we used N = 4.5 x 
10 3 points of the time sequence of p m and p m , averaged 
over 5 x 10 4 ICs. At small H e ff = 0.0001 we observe 
in the parameter space for T>{p : p) that the structures 
Hi, I? 2 , C_ i, F can be clearly recognized (they have sharp 
borders). For parameters inside the ISSs the quantum- 
classical dynamics is strongly correlated, while outside 
the ISSs the correlation goes to zero. Inside the ISS Bi we 
observe that correlations also go to zero, but its border is 
still well defined. For H e ff = 0.082 the result is inverted, 
as inside most part of the ISSs the correlation is zero and 
outside it remains almost the same as from Fig. [3 (a). 
In this case the ISS B 2 disappeared, the C_i and F yet 
exist with sharp borders, and some remaining borders of 
B\ are visible. These results demonstrate that the sharp 
borders of the ISSs tend to survive in the parameter space 
(for reasonable h e ff) when V(p,p) is plotted, while the 
ISSs are already blurred and unrecognizable when the 
ratchet current is plotted [compare Fig. |3jb) to Figs. [3 a) 
and (b)[. In other words, the way the ISSs disappear in 
both cases is fundamentally distinct. It is obvious that 
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Figure 3. (Color online) Plotted is V(p,p) for (a) h e ff = 
0.0001 and (b) h e ff = 0.082, and the activation/suppression 
8{h ef f) for (c) h e ff = 0.01 and (d) h e ff = 0.082. 

using huge values for H e ff the correlation in the whole 
parameter space will go to zero. 


C. Quantum activation of the ratchet current 

Another interesting related phenomenon is that the 
quantum fluctuations not only lower the Cs, but also ac¬ 
tivate them. In order to show this we calculate 5(h e ff) = 
|Cs(?ie//)| — |Cc(Te// = 0 )| and plot a map of the cur¬ 
rent activation/suppression in the parameter space. For 
5 > 0 activation of the Cs occurs while for 5 < 0 the Cs 
decreases due to vacuum fluctuations. Results are shown 
in Fig. [3])c) and (d) for the same parameter space from 
Fig. 0 We use green to dark-green (gray to dark-gray) 
colors for increasing activations of the Cs and white, red, 
blue to black (white, dark-gray to black) for increasing 
suppressions of Cs- Surprisingly there are many regions 
in the parameter space where activations occur. The 
physical origin of the activations depend on the specific 
dynamics at a given parameter combination. Along the 
diagonal dark-green line (below the ISS B i) where large 
activations (A ~ 6 ) occur in Fig. [3jc) , we conclude from 
classical results m that the activation is due to a crisis 
bifurcation. Additional numerical results revealed that 
the activation close to K = 2.5 and 7 = 0.85 is related 
to a symmetry breaking of the multistable attractors. 


IV. CONCLUSIONS 

Concluding, the quantum ratchet current was ana¬ 
lyzed in the parameter space of the asymmetrically kicked 
rotor. We have shown that isoperiodic stable struc¬ 
tures with sharp borders play an essential role in the 
quantum-classical transition. Quantum currents were ob¬ 
tained numerically exact by solving the Markov stochas¬ 


tic Schrodinger Eq. (|2j) (Cq) , and in a semiclassical ap¬ 
proach from the semiclassical map (j3j(Cs). For the clas¬ 
sical currents (Cc) we used the map which was de¬ 
rived via direct integration of the kicked Langevin equa¬ 
tion with a white noise. The quantum-classical transition 
was studied by means of a comparison between the Cs 
and Cc- We show how and to which extent the classi¬ 
cal generic isoperiodic stable structures mmm are 
destroyed by quantum fluctuations. Even though these 
structures become blurred and start to disappear in the 
quantum world, their sharp borders tend to persist when 
the correlation between quantum and classical currents is 
investigated. For this aim we have used the distance cor¬ 
relation measure [22]. Remarkably we have found a sta¬ 
tistical dependence of the quantum and classical currents 
when parameters are chosen at the borders of the isope¬ 
riodic stable structures. Moreover, very often this de¬ 
pendence extends to the interior of the isoperiodic struc¬ 
tures. In addition, we calculated the quantum activa¬ 
tion/suppression of the current in the parameter space. 
Surprisingly many regions were found where quantum 
fluctuations enhance the Cc- Even though the distance 
correlation and the current activation are not directly 
observed for one simulation or experiment, they can be 
obtained by comparing two simulations in distinct semi¬ 
classical regimes and allow us to understand better the 
quantum-classical transition of ratchet currents. 

It is generally accepted that dissipative quantum dy¬ 
namics inducing decoherence will blur any quantumness 
and eventually will lead to a dynamics that may well be 
described in classical terms. In this paper we show, how¬ 
ever, that for the ratchet current in a zero-temperature 
environment, this quantum-classical transition is far from 
trivial and depends crucially on the chosen parameters. 
We prove that a replacement of quantum fluctuations 
by classical thermal fluctuations is overly simplistic and 
cannot be used for the whole parameter space uniformly. 
It is clear that a single quantity like the Cq can be 
matched by a Cc if one allows an additional free pa¬ 
rameter (here temperature) to vary. Crucially, however, 
different points in parameter space require very differ¬ 
ent temperatures for this fit procedure (see Fig. 2), and 
in this sense a quantum-classical correspondence, where 
quantum fluctuations are replaced by thermal fluctua¬ 
tions, cannot exist. In addition, there is no one scaling 
which gives the quantum-classical correspondence in the 
whole parameter space. One of the reasons for this be¬ 
havior can be traced back to the complicated dynamics 
in the case of mixed phase space structure. In our con¬ 
text we can write Cs = y~); cqC^, where the sum is 
over the attractors, Cg 1 ^ is the ratchet current from at¬ 
tractor i weighted by the area of the corresponding basin 
of attraction, and 07 = /x,Tj is the statistical weight of 
each attractor, where 77 is the mean lifetime on attrac¬ 
tor i and /x,; = izm n _ > . 00 (mi)/n where rrii is the number 
of times a given trajectory visited the attractor i during 
the time n. Quantum fluctuations allow the trajectory to 
jump between attractors and 07 is the quantity which de- 
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pends, in a nontrivial way, on h e ff and may also depend 
on the size/shape of the basin of attraction. The degree 
of quantum-classical correspondence is therefore strongly 
dependent on the local dynamics, and general statements 
about the scaling of a* are difficult. However, our results 
indicate a qualitatively different correspondence for pa¬ 
rameters inside and outside the ISS regions. Indeed, we 
have shown that the way the ratchet current responds to 
the quantum fluctuations depends sensitively on whether 
the parameters are chosen within the ISSs’ limits. It indi¬ 
cates that the classical-quantum transition is much more 
abrupt inside the ISSs, at least regarding the ratchet cur¬ 
rent. 

Our results are an important step towards the goal 
to experimentally observe quantum isoperiodic struc¬ 
tures in the context of ratchet effects with cold atoms 
[2B1 27]. Moreover, we strongly believe that quantum ac¬ 
tivations deserve further analysis, both theoretically and 


experimentally. In general, the analysis of the complex 
quantum-classical transition in the parameter space (in¬ 
cluding the isoperiodic structures), is not restricted to 
ratchet currents. Other observable of a quantum system, 
whose classical counterpart presents regular and chaotic 
dynamics, could be used. We mention for example the 
motion of atoms in an optical lattice |2Sl ]2U] and the 
electronic wave-packet dynamics in Rydberg atoms sub¬ 
jected to a strong static magnetic field j2TjjH32!- 
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